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Abstract. Experiments with gold-gold collisions at RHIC have revealed (i) 
stronger suppression of charmonium production at forward rapidity than at 
midrapidity and (ii) the similarity between the suppression degrees at RHIC and 
SPS energies. To describe these findings we employ the model that includes 
nuclear shadowing eff'ects, calculated within the Glauber-Gribov theory, rapidity- 
dependent absorptive mechanism, caused by energy-momentum conservation, and 
dissociation and recombination of the charmonium due to interaction with co- 
moving matter. The free parameters of the model are tuned and fixed by 
comparison with experimental data at lower energies. A good agreement with 
the RHIC results concerning the rapidity and centrality distributions is obtained 
for both heavy Au+Au and light Cu+Cu colliding system. For pA and A+A 
collisions at LHC the model predicts stronger suppression of the charmonium and 
bottomonium yields in stark contrast to thermal model predictions. 



1. Introduction 

The investigation of nuclear matter under extreme temperatures and densities, and 
the search for a predicted transition to a deconfined phase of quarks and gluons, the so- 
called Quark-Gluon Plasma (QGP), is one of the main goals of heavy-ion experiments 
at ultrarelativistic energies. Both theorists and experimentalists are looking for 
genuine QGP fingerprints, that cannot be masked or washed out by processes on 
a hadronic level. Charmonium was proposed about two decades ago [1] as one of the 
most promising QGP messengers because its yield would be significantly suppressed 
due to color Debye screening in the plasma phase. Also, due to the small interaction 
cross section of J/ijj in hadronic matter, charmonium spectrum is expected to carry 
information about the early hot and dense stage of nuclear collision. Since the volume 
of the produced QGP depends on the collision energy, centrality and mass of colliding 
nuclei, it is generally believed that the suppression of the J/V' yield would increase 
with rise of the aforementioned factors. 

Therefore, the PHENIX measurement of the nuclear modification factor of 
charmonium at top RHIC energy = 200AGcV [2] uncovered at least two 
unexpected features. Firstly, compared to charmonium suppression in lead- lead 
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collisions at SPS energy Eiab — 160AGeV [3] the level of suppression at midrapidity 
at RHIC was found to be quite similar for the same number of participants despite 
the order of magnitude difference in center-of-mass energies of heavy-ion collisions. 
Secondly, J/t/i suppression increases unambiguously with rising rapidity, whereas the 
highest energy density and the most dense medium should be produced at y = 0. 
These findings, as well as the experimental results for dAu collisions, where the plasma 
formation is very unlikely, attract the attention to the whole variety of the processes, 
including both initial and final state effects, that are responsible for the charmonium 
production and its propagation through hot and dense medium. 

The paper is organized as follows. Section [2] presents a description of the model 
that contains a comprehensive treatment of initial-state nuclear effects, such as nuclear 
shadowing and nuclear absorption, and final state interactions with the co-moving 
matter. Comparison with the available experimental data and predictions for Pb + Pb 
collisions at LHC energy is given in Sec. [3] Finally, conclusions are drawn in Sec. 31 



2. Description of the model 

Nuclear effects in nucleus-nucleus collisions are usually expressed through the so-called 
nuclear modification factor, i?^'^'(6), defined as the ratio of the J/ip yield in A+A and 
pp scaled by the number of binary nucleon-nucleon collisions, n(b). We have then 
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Here (7AB{b) = 1 — exp[—app AB TAB{b)], TAB{b) = / d?' sTA{s)TB{b — s) is the nuclear 
overlapping function, TA{b) is obtained from Woods-Saxon nuclear densities, and 

= <jppABTA{s)TB{b- s)/aAB{b) . (2) 

where upon integration over d?s we obtain the number of binary nucleon-nucleon 
collisions at impact parameter 6, n{b). The three additional factors in the numerator 
of Eq. IH), 5'*'', 5'"'"' and S*™, denote the effects of shadowing, nuclear absorption, and 
interaction with the co-moving matter, respectively. Let us discuss them briefly. 

The nuclear absorption is usually interpreted as suppression of J/ -0 yield because 
of multiple scattering of a cc pair within the nuclear medium. At low energies the 
primordial spectrum of particles created in scattering off a nucleus is mainly altered by 
(i) interactions with the nuclear matter they traverse on the way out to the detector 
and (ii) energy-momentum conservation. For A+A collisions these effects can be 
combined into the generalized suppression factor 

[l-cM-i{x+)^QQATA{b))] 



i{x+)i{x-)al^ABTA{s)TB{b - s) 

[l-eM-i{^-)<^QQBTB{b-s))] , (3) 



where x± = (^a;| — 4M2/s ± x-p)/'2., and £,{x±) = (1 — e) -I- exj determines both 
absorption and energy- momentum conservation. In [4] it has been found that 7 = 2, 
e — 0.75 and ctqq — 20 mb give a good description of data. This corresponds to 
Cabs = 5 mb at mid-rapidity and agrees well with other studies. 
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Secondly, coherence effects will lead to nuclear shadowing for both soft and hard 
processes at RHIC, and therefore for the production of heavy flavor. Shadowing can 
be calculated within the Glauber-Gribov theory [5] , and we will utilize the generalized 
Schwimmer model of multiple scattering [6] . In this case the second suppression factor 
in Eq. ((T|) reads 

^''^^'"'^^ " 1 + AF{yA)TA{s) 1 + BF{yB)TB{b~s) ' 
where the main contribution to the function F{y), that encodes the dynamics of 
shadowing, comes from the gluon rather than from the quark shadowing [7,8]. At 
SPS energy the nuclear absorption dominates over the shadowing, whereas RHIC 
already belongs to the high-energy regime. Nuclear shadowing is non-negligible at mid- 
rapidity, and the combined effect of shadowing and energy-momentum conservation 
should be accounted for at forward rapidities. At LHC shadowing will be very strong 
even at y = 0, while energy-momentum conservation becomes a minor effect. 

Finally, the processes of dissociation and recombination of cc pairs in the dense 
medium should be taken into account. We employ the co-movers interaction model 
(CIM) [9] that was recently modified to incorporate the recombination mechanism 
into consideration [10]. Assuming a pure longitudinal expansion and boost invariance 
of the system, the rate equation which includes both dissociation and recombination 
effects for the density of charmonium at a given production point at impact parameter 
s reads 

dNj/^{b,s,y) 
^ Tr 

- Nc{h,s,y)N-,{h,s,y)\ , 

where A^™, Nj/^, and A"c(c) is the density of comovers, J/ip and open charm, 
respectively, and CTco is the interaction cross section for both dissociation of 
charmonium with co-movers and regeneration of J from cc pairs in the system 
averaged over the momentum distribution of the participants. It is the constant of 
proportionality for both the dissociation and recombination terms due to detailed 
balance Nj/^{b, s,y) N'^°{b, s,y) = Nc{b,s,y) Nc{b,s,y). The solution of Eq. ^ can 
be approximated by 

S^''{b,s,y)= exp{-a,o[A^^°(6,s,2/) 



N-''ib,s,y)Nj/^(b,s,y) (5) 
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where 



N' 
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C{y) = ' . (6) 

Details of the model can be found in [10]. The quantities in Eq. ([6]) are all related to 
pp collisions at the corresponding energy and are taken from experiment. Note, that 
the extension of the CIM to the recombination effects does not imply any additional 
parameters. 



3. Heavy quarkonium at RHIC and LHC 



The density of open and hidden charm at mid-rapidity in pp collisions at ^/s = 200 
GeV has been reported in [11]. In the left picture of Fig. [T] we present the results of 
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Figure 1. Results for J/ip suppression in Au+Au at RHIC {y/s = 200 GeV) at 
mid- (left figure), and at forward rapidities (right figure). Data are from [2], The 
solid curves are the final results. The dash-dotted ones are the results without 
recombination (C = 0). The dashed line is the total initial-state effect. The 
dotted line in the right figure is the result of shadowing. In the left figure the last 
two lines coincide. 



our model compared to experimental data at mid-rapidity. The different contributions 
to J/ip suppression arc shown. Note that at mid-rapidities the initial-state effect is 
just the shadowing. As discussed above nuclear absorption due to energy-momentum 
conservation is present at forward rapidities but is negligibly small at mid-rapidities. 
One can see that the J/4> suppression at forward rapidity is somewhat larger than 
that at mid-rapidities, in full accord with experimental data. This is due to both the 
recombination term and the initial-state effects. The latter are stronger for forward 
rapidities. 

For consistency, we have also made calculations for the J/ip suppression in 
Cu+Cu collisions at RHIC using the same parameters as above for Au+Au collisions. 
The results are shown in Fig. [21 and are in good agreement with the experimental 
data, except maybe for peripheral collisions, where the error bars are quite large. 
Concluding, our procedure gives a reasonable description of data both at mid- and 
forward rapidity for different collision systems at RHIC. The effect of recombination 
is more pronounced at midrapidity. 

Based on our previous discussion, it is obvious that dissociation-recombination 
effects will be of crucial importance in Pb+Pb collisions at LHC {^/s = 5.5 TeV). 
Assuming that the energy dependence of open charm and J/^ in pp collisions is the 
same (between RHIC and LHC energies), the energy dependence of the parameter 
C will be that of app/app. The total and differential cross section for charm can be 
calculated using perturbative techniques [13]. The calculations for low energies are in 
agreement with data, yet predictions for RHIC and Tevatron energies are lower than 
the data. Therefore, the extrapolation to LHC is quite uncertain. If we parameterize 
the energy dependence of open charm production as cr'^'^ oc s", with a = 0.3 and use 
the values of non-diffractive app as 34 mb for RHIC and 59 mb for LHC, we obtain 
C = 2.5 at LHC - a value about four times larger than the corresponding one at 
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Figure 2. Results for J/ij) suppression in Cu+Cu at RHIC {y/s = 200 GeV), at 
mid- (left) and at forward rapidities (right). For details, see caption of Fig. [T] 
Data are from [12]. 



RHIC. In view of that we consider that reahstic values of C at LHC are of the range 
2 to 3. In Fig. [Sja) we have calculated the J/ip suppression at LHC for several values 
of C, including the case of absence of recombination effects (C ~ 0). Although the 
density of charm grows substantially from RHIC to LHC, the combined effect of initial- 
state shadowing and comovers dissociation appears to overcome the effect of parton 
recombination. This is in sharp contrast with the findings of [14], where a strong 
enhancement of the J/'ijj yield with increasing centrality was predicted. Note that in 
our approach only comovers and open charm produced at the same impact parameter 
as the initial J/tp are allowed to interact, whereas models assuming global equilibrium 
of the produced charm with the medium allow for recombination of cc pairs from the 
whole volume of the fireball. 

Finally, we would like to discuss the impact of initial state effects on bottomonium 
production. The absorptive cross section for T is 40-50% smaller than that for J/ij; 
and "0', and energy- momentum conservation mechanisms are pushed to higher xp 
due to the large mass of the bottomonium. Therefore, nuclear absorption for T at 
LHC is expected to be quite small. The suppression of bottomonium due to gluon 
shadowing in Pb-|-Pb collisions at LHC is shown in Fig. ^h) for several rapidities. 
The suppression is about 50% from mid-central to central collisions, and would be 
the same for all members of the T family. This establishes the baseline for further 
calculations of bottomonium dissociation and recombination in the final state. 

4. Conclusions 

The effects of recombination of cc pairs into J/ip are incorporated in the comovers 
interaction model. These effects are negligible at low energies (SPS) due to the low 
density of open charm. The model does not assume thermal equilibrium of the matter 
produced in the collision and includes a comprehensive treatment of initial-state 
effects, such as shadowing, nuclear absorption and energy-momentum conservation. 
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Figure 3. Left: J/i/' suppression in Pb+Pb at LHC {\/s = 5.5 TeV) at mid- 
rapidities for different values of the parameter C. The upper Hne is the suppression 
due to initial-state effects (shadowing). Right: Centrality dependence of T 
suppression due to gluon shadowing in Pb+Pb collisions at = 5.5 TeV. 



In our approach, the magnitude of the recombination effect is controUed by the 
total charm cross section in pp colhsions. Using it as an input, the centrahty and 
rapidity dependence of experimental data is reproduced both for Au+Au and Cu+Cu 
collisions at full RHIC energy. For LHC we are lacking the experimental information 
and should rely on estimates from theoretical models. For a reasonable choice of 
parameters, we predict that the suppression observed at RHIC and lower energies 
will still dominate over the recombination effects. This is due to the large density of 
comovers and to the strong initial-state suppression at these ultra-relativistic energies. 
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